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The capability to probe the dispersion of elementary spin, charge, orbital, and lattice excitations 
has positioned Resonant Inelastic X-ray Scattering (RIXS) at the forefront of photon science. Here 
we develop the scattering theory for RIXS on superconductors, calculating its momentum-dependent 
scattering amplitude. Considering superconductors with different pairing symmetries we show that 
the low-energy scattering is strongly affected by the superconducting gap and coherence factors. 
This establishes RIXS as a tool to disentangle pairing symmetries and to probe the elementary 
excitations of unconventional superconductors. 
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Introduction In the past decade, resonant inelastic 
x-ray scattering (RIXS) [l|, |2| has made remarkable 
progress as a spectroscopic technique, establishing it- 
self as an experimental probe of elementary spin [3|-l8j, 
orbital 0, [To| , and lattice excitations [ll|. In quite a 
number of cases, theoretical considerations have preceded 
and stimulated these experimental advances, prominent 
examples being the theoretical demonstration of the 
presence of strong single-magnon scattering channels in 



cuprates [12|, [13[ and iridates [Ijj. Being a photon- 
in/photon-out spectroscopy, both the energy hu and the 
momentum change q of the scattered photon are mea- 
sured. As the energy and momentum lost by the photon 
are transferred to intrinsic excitations of the material un- 
der study, direct information on the dispersion of those 
excitations becomes available. The resonant character of 
the technique is due to the energy of the incident photon 
being chosen such that it coincides, and hence resonates, 
with an x-ray absorption edge of the system [l], Q . This 
year the energy resolution of RIXS has reached ~30 meV 
in the hard x-ray regime [7j , will reach 50 meV in the soft 
x-ray regime by building on present instrumentation (l5| 
and is designed to reach 11 meV at the Cu L-edges at the 
NSLS-II presently under construction [16]. This brings 
the RIXS energy resolution well into the regime of the en- 
ergy gap of cuprate superconductors, which stretches out 
to 119 meV for mercury-based high T c systems [Tjj. Con- 
sequently the fundamental question arises of how the su- 
perconducting (SC) state leaves its fingerprints in RIXS 
spectra — in particular whether and how RIXS is sen- 
sitive to the phase and amplitude of the SC gap and to 
quasiparticle excitations. 

Probing the order parameter in unconventional super- 
conductors is generally the first step for an investigation 
of the pairing mechanism and of the character of the SC 
state. Compared to the available spectroscopic methods, 
such as scanning-tunneling spectroscopy (STS), photo- 
emission spectroscopy, optical spectroscopy or inelastic 
neutron scattering, RIXS uniquely combines the advan- 



tages of bulk-sensitivity and availability of momentum 
resolution while at the same time requiring only small 
sample volumes. Here we show how the sensitivity of the 
RIXS process to the dynamical structure factor (DSF) 
of the electron spin and density in the SC state enables 
the investigation of SC quasiparticle excitations. In par- 
ticular, we find that the momenum dependence of RIXS 
spectra is intrinsically determined by the pairing symme- 
try, being sensitive not only to the magnitude of the SC 
gap and to the presence of nodes on the Fermi surface 
but also to the phase of the order parameter. This phase 
sensitivity is due to the appearance of SC coherence fac- 
tors which, for instance, in STS determine to large extent 
the quasiparticle interference in the presence of impuri- 
ties [li-lii. 

Dynamical structure factors in RIXS In what follows 
we concentrate on the so-called direct RIXS process at 
the transition metal (TM) ion L edges, in which the in- 
coming photon resonantly excites the core shell 2p elec- 
tron into the 3d shell which consequently decays into an 
outgoing photon and a charge, spin, or orbital excitation 
in the electronic system [2|. In this case the RIXS cross 
section can be written as 




7 e (q,^) = ^|(/|O q , € 

f 



i)\ 2 5(hcu-E f + E z ), (1) 



where \i) (|/)) is the initial (final) state of RIXS pro- 
cess with energy Ei (Ef) and hu (q) is the transferred 
photon energy (momentum). Furthermore, in the fast 
collision approximation (FC A) 0, EH, 0, Hsj the Fourier 
transformed RIXS transition operator O q , e can be writ- 
ten as O q5 e = W^Pa + Wg • S q , where the operators 

p q = Xlkcr ^k+qcr^kcr an d are the density and spin of 
conduction electrons j25l427| . The so-called RIXS form 
factors W£ and depend on the TM ion, the specific 
geometry of the experiment, and on the polarization e 
of the incoming and outgoing photon — their precise de- 



pendencies are provided in Refs. LL3|, [25|, [26|. In terms of 
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these form factors the RIXS cross section is 



(a) anisotropic s wave 
-7T 



(c)q vector 



j e (q,w) = |w e c rx c (q,w) + |wj|V(q,w), (2) 

where x c = J2f l(/l/°qK)| 2< K^ ~ Ef + ^) * s the charge 
and x s = £/ |(/|5*|z>| 2 *(fia;- J57/ + the spin dynam- 
ical structure factor (DSF), where one assumes that the 
spin DSF has the same momentum and energy depen- 
dence for all the three different components of the spin 
operator — as is the case of unconventional supercon- 
ductors (3l|. As is clear from the above, the amplitude 
of the RIXS form factors W£ and can be tuned by 
properly adjusting the experimental conditions in RIXS. 
This implies that RIXS at L-edges can measure either 
spin or charge DSF depending on the chosen polariza- 
tion, which is a unique feature of RIXS spectroscopy. 
Note that the FCA has been successfully used to describe 
low energy excitations in RIXS, e.g., at the Cu L edge 
of both undoped and doped cuprates (cf. theoretical cal- 
culations [HI, HH, and experimental results with their 
interpretation fully based on the FCA H, [To|, [HI), at 
the Fe L edge of the iron arsenides [29|, and at the Ir 
L edge in iridates [30]. It has also been shown theoreti- 
cally that when the incoming photon energy in RIXS is 
tuned to a TM ion resonant edge in a TM oxide, the FCA 
describes the RIXS spectrum well (25| . 

In the following we concentrate on determining the 
properties of DSF for different types of singlet-pairing 
superconductors, and to be even more specific below we 
consider the case of Cu ions in lattices with tetragonal 
symmetry, i.e., the one which has direct relevance to the 
high T c superconductors. Our main aim in this context is 
to establish how a variation of the phase of the SC order 
parameter is reflected in the spin and charge DSF. Fol- 
lowing the most direct theoretical inroad and avoiding 
model-specific technical details, we performed calcula- 
tions for a singlet-pairing superconductor with a SC order 
parameter varying along the Fermi surface. Even if by 
this electron correlations are not fully taken into account, 
this approach is commonly used - and is very successful 
to calculate quasiparticle interference in cuprates [2ol — 



22|, |32|. Besides this, in the supplemental material we 
show that it is actually possible to introduce strong cor- 
relations between electrons into the calculations and that 
such does not affect the main results on the SC electronic 
system presented below. 

DSF for a superconductor Quasiparticle excitations 
in a single band superconductor with singlet-pairing are 
described by the Hamiltonian H—£pN = E^ 7^7^ > 
where = \j s\ + l^k| 2 is the quasiparticle energy 
dispersion depending on the SC gap function Ak and 
on the dispersion of the bare electrons e^. The Bo- 
goliubov quasiparticle operators 7 kcr are related to the 
d electron operators via d^ = w k 7kt — ^k7lk| anc ^ 

d U = <7k| + ^k7l k t> with K| 2 Ok| 2 ) = \ (1 ±e k /£k) 
and u\v\^ = ^A^/E^. The Bogoliubov transformation 




FIG. 1: (color online) Order parameters over the Brillouin 
zone for an anisotropic s wave (a) and a d wave (b) supercon- 
ductor with an isotropic Fermi surface (solid line) . The nodal 
lines (dashed) of the order parameter (i.e., where Ak = 0) 
intercept the Fermi surface in the nodal points (kjv). On the 
other hand, anti- nodal points (k^) are defined by the points 
on the Fermi surface where the order parameter has maxima. 
In the d wave case one can distinguish between two kind of 
particle- hole excitations (c), i.e., the ones with the same sign 
of the order parameter, and the one with sign reversal, with 
a transferred momentum connecting two regions of the Bril- 
louin zone with a different sign of the order parameter. 



allows the evaluation of the DSF for a SC by evaluat- 
ing the transition amplitudes (/|p q |i) and (f\S^\i) be- 
tween the ground state \i) and any excited state |/) of 
the Hamiltonian. At zero temperature, the excited states 
which contribute to DSF have the form 7k+ q a ^-h -o\9) 
with a transition energy of £?k+q + ^-k- Using the Bo- 
goliubov transformation one then finds that the DSF for 
a superconductor reads 



X c ' s (q )W ) = £ 



1 ± 



Re ( A k A k+q) T £k£k+q 



-Ek-Ek+q 
X 5(hLJ - -Ek+q - Ek) 



(3) 



where ± sign is for the charge (spin) DSF [33-35]. Thus, 
the DSF is a sum over all momenta within the Brillouin 
zone, where the transition amplitudes are strongly influ- 
enced by the character of the SC state. Although quasi- 
particle interactions substantially affect the DSF, they do 
not alter its intrinsic sensitivity to the character of the SC 
state (in particular to the symmetry of its gap function), 
as can be seen, e.g., at the random-phase approximation 
level [36], or by considering a strongly correlated system 
with Hubbard interactions, see supplemental material. 

To determine in detail how the RIXS spectra of un- 
conventional superconductors reflect the phase of the or- 
der parameter, we consider (i) a d-wave pairing with 
Ak oc (cosk x — cos k y ) and (ii) an anisotropic s- wave 
pairing with Ak oc | cos k x — cos k y | , i.e., two pairing sym- 
metries which differ from each other only in the SC order 
parameter phase. Maps of the considered gap functions 
and the Fermi surface are shown in Fig. [TJ In the s wave 
case, the DSF x c (q^) 7 due to the sign in Eq. ((3]), 
is non-zero all over the Brillouin zone, while in the d 
wave case excitations combining momenta with opposite 



3 



(a)isotropic FS (b)nested FS 




012345012345 



FIG. 2: (color online) Dynamical structure factor (DSF) 
in a SC for an isotropic (a) and for a nested Fermi sur- 
face (b) with transferred momentum q = (/cf, /cf)/v / 2 and 
q = (7r, 7r) respectively, for charge excitations (% c ) with 
anisotropic s wave (solid line) and d wave (dashed line), or 
alternatively for spin excitations (x s ) with d wave (solid line) 
and anisotropic s wave (dashed line) order parameter. In the 
d wave (anisotropic s wave) case, charge (spin) DSF of sign 
reversing (sign preserving) excitations vanish and the spec- 
tral weight is strongly suppressed, while in the anisotropic s 
wave (d wave) case no phase cancellation occurs, allowing the 
coherence peak (b) to appear in the case of a perfectly nested 
Fermi surface. 



phases of the order parameter, i.e., sign reversing pro- 
cesses, are suppressed. Note that for the spin DSF this 
situation is inverted. 

Phase sensitivity The calculated quasiparticle spec- 
tra for the above two pairing symmetries are shown in 
Fig.[2fa) for a fixed and relatively large momentum trans- 
fer q — fc_p, where kp is the Fermi radius. A detailed ex- 
pansion of the DSF for q <C kp is provided in the supple- 
mentary materials. Note that direct RIXS at the L-edge 
in 2D cuprates allows momentum transfers q < 0.877T [28| 
and therefore one is able to access momentum transfers q 
of the order of /cp- Such momentum transfers correspond 
to a momentum vectors combining two anti-nodal points 
on the Fermi surface [see inset of Fig. Efa)] with the 
same gap value Ak A = A. For such an excitation, the 
sign of the order parameter in the s wave case is pre- 
served whereas in the d wave case it is reversed. Clearly 
the spectral weight at hu = 2|A| in the charge DSF is 
enhanced in the s wave case. On the other hand, for 
the d wave case the sign reversal leads immediately to a 
suppression of the DSF according to Eq. (j3]). Therefore, 
the spectral weight at Huj = 2|A| is expected to be very 
small which is confirmed by the dotted curve in Fig.[2fa). 
On the other hand, spin DSF is suppressed in the s wave 
case, while the sign reversal enhances the spectral weight 
at Hlj = 2|A|. Note that due to the gap being equal in 
magnitude for both cases, the obtained effect is entirely 
due to phase changes of the SC order parameter along 
the Fermi-surface. 

In unconventional superconductors, where the pairing 
is generally considered as mediated by antiferromagnetic 
spin fluctuations [37], the SC gap function is expected 
to exhibit a sign reversal between the Fermi momenta 



connected by characteristic wave vector Q of the spin 
fluctuations [38|. The conduction electrons of such su- 
perconductors show a tendency to Fermi-surface nesting 
with a typical nesting vector of Q = (n, tt). In Fig. Gfb) 
the scattering intensity as function of energy for the two 
pairing symmetries considered above are shown again, 
but now for a perfectly nested cuprate-like Fermi surface 
(see inset), with a transferred momentum equal to the 
nesting vector. It is clear that, in comparison to the case 
of an isotropic Fermi surface, the coherence peak in the 
anisotropic s wave pairing case appears now strongly en- 
hanced due to the nesting effect. As in the case of an 
isotropic Fermi surface, sign reversing excitations occur- 
ring in the d wave case in the charge DSF are strongly 
suppressed, as well as sign preserving excitations in the 
anisotropic s wave case in the spin DSF. 

To highlight its strong dependence on the order pa- 
rameter phase, we show in Fig. [3] the DSF for a fixed 
energy Hlj = 2|A| as a function of momentum q in the 
entire Brillouin zone, both for the anisotropic s wave 
and the d wave pairing, for an isotropic and a perfectly 
nested Fermi surface. For an isotropic Fermi surface, 
d wave (anisotropic s wave) pairing strongly suppresses 
the spectral weight in the charge (spin) DSF for trans- 
ferred momenta that correspond to sign reversing (sign 
preserving) excitations, while no suppression occurs in 
the anisotropic s wave (d wave) case (Fig. [3j first row). 

For a nested Fermi surface, coherence peaks are clearly 
visible in the charge (spin) DSF in the anisotropic s 
wave (d wave) case for any of the nesting vector q = 
(±7r, ±7r), while they are strongly suppressed in the d 
wave (anisotropic s wave) case (Fig. 02 second row). 
Clearly the symmetry of the order parameter is reflected 
by the symmetry of the charge and spin DSF spectrum. 
Since the charge and spin DSF are complementary re- 
spect to the spectral suppression of the sign reversing 
and sign preserving excitations, the phase sensitivity is 
enhanced when these two components are fully disen- 
tangled. This can be done by tuning the polarization 
dependence in the form factors W£ and in Eq. ((2]). 

Origin of the phase sensitivity It occurs that a strong 
dependence on the SC order parameter phase in DSF is 
found for transferred momenta q > 0.1 kp and when the 
transferred energy is close to twice the energy of the SC 
gap Hlj 2|A|. This sensitivity to the SC order parame- 
ter phase can be better understood if we further confine 
our discussion to the case where not only the transferred 
momentum is rather large, i.e., q ~ kp, but also the 
momentum k is on the Fermi surface. The main con- 
tributions to the DSF correspond in fact to those exci- 
tations close to the Fermi surface which fulfill the con- 
ditions £k <C Ak and £k+ q <C Ak+ q . Assuming an un- 
conventional superconductor with a pairing governed by 
a phase dependent order parameter Ak = |Ak|e*^ k the 
DSF in Eq. (j3]) for excitations near the Fermi surface 
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FIG. 3: (color online) Charge (% c ) and spin (% s ) DSF at fixed energy ftw = 2|A|, as a function of transferred momentum 
(first two and last two columns respectively) for an anisotropic s wave (first and fourth columns) and a d wave (second and 
fifth columns) SC, with an isotropic (first row) and a nested Fermi surface (second row). In the middle column, the total DSF 
is shown, given by the sum of charge and spin contributions (x° + X s ) 5 f° r both an anisotropic s wave and a d wave SC, with 
isotropic (first row) and a nested Fermi surface (second row). The order parameter symmetry shows up in the momentum 
dependence of both charge and spin DSF, in particular for a perfectly nested Fermi surface, where phase cancellations suppress 
the coherence peak at q = (±7r, ±7r) in the charge DSF in the d wave pairing case, and in the spin DSF in the anisotropic s 
wave case. The phase sensitivity is lost instead when one sums up the two contributions (x° + X a ) 5 with no differences respect 
to the order parameter phase. The isotropic and the nested Fermi surface spectra are plotted with different intensity scales. 



(Jhjuj < 2 A) becomes 

X c ' s (q, M « I 1 ± cos (^k - 0k+q)] 

x<K^-|A k+q |-|A k |). (4) 

Thus, the momentum-dependent intensity distribution of 
the low energy DSF mainly represents the variation of the 
SC order parameter phase along the Fermi surface. 

Comparison with other spectroscopies In principle 
also other two-particle spectroscopies (see, e.g., Ref. [2| 
for an overview) can be directly sensitive to the DSF of 
superconductors. Even if none can match RIXS in mea- 
suring both spin and charge DSF of superconductors, al- 
ready probing either of the two is in general challeng- 
ing as in, e.g., electron energy loss spectroscopy (EELS) 
one cannot reliably measure spectra with high momen- 
tum transfers, inelastic neutron scattering (INS) does not 
directly probe the charge DSF and non-resonant inelas- 
tic x-ray scattering (NIXS) is extremely photon- hungry. 
Nevertheless, transition amplitudes of the same type as 
in Eqs. (j3|4p are also encountered when determining the 
scattering rate of conduction electrons in the presence of 
impurities, as observed in the surface-sensitive STS [lil - 
This is because Eqs. (j3|4p have a similar structure 
as the ones which are known to govern the quasiparti- 
cle interference [in which case the transition amplitudes, 
whose sum over the momentum k contribute to the DSF 
in Eq. (|3j), are termed 'coherence factors'] in the pres- 
ence of impurities. Since the quasiparticle interference 



patterns explored by STS have turned out to be very suc- 
cessful in uncover the pairing symmetries of the uncon- 
ventional SC [l8l-[23|, this gauges the potential of RIXS 
to observe and unravel symmetries of SC pairing and 
pairing- mediators. Compared to STS, however, RIXS 
has a succinct conceptual advantage. Whereas the theo- 
retical interpretation of STS in the framework of quasi- 
particle interference relies crucially on the form of the 
underlying impurity system showing various components 
of scattering [24] , in the case of RIXS the interpretation 
of spectroscopic features does neither rely on the presence 
of impurities in the superconductor nor on the modeling 
thereof. 

Conclusions We have shown that RIXS, in contrast 
to other well-known two-particle spectroscopies, is di- 
rectly sensitive to the spin and to the charge dynamical 
structure factor (DSF) of a superconductor. In partic- 
ular we have shown that the DSF of a superconductor 
observed in RIXS is very sensitive to the symmetry of 
the order parameter. This is rooted in the quasiparti- 
cle spectra reflecting sign-reversing excitations at large 
transferred momenta which arise for order parameters 
with a phase that varies over the Fermi surface. This, 
together with the recent experimental successes of RIXS, 
including in particular the major enhancements in resolu- 
tion and pioneering study of hole doped cuprates [6| , es- 
tablishes the potential of RIXS as a versatile and practi- 
cal spectroscopic technique to investigate the fundamen- 
tal properties of superconducting materials. 
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SUPPLEMENTAL MATERIAL 



transition energy i^k A +q + is approximately 



Dynamical structure factor for small (q <C /cf) 
scattering 



5 
4 

< 3 

3 

«2 2 
1 




(a)s wave 



(b)d wave 



0.02*^(1,1) 



(0,0) 0.02* F (1,0) 0.02it f (l,l) 



(0,0) 0.02^(1,0) 



£ kA+q + £ kA ^/|A|2+(^^) . (5) 



where Vk = de^/hdk is the electron group velocity. Thus, 
as it is seen in Fig.H^b) for E > 2 A we obtain a coherence 
peak showing the known dispersion of Fig.H^a), but with 
less spectral weight since the excitations are restricted to 
the vicinity of the anti-nodal points. The q-dependence 
of the excitation energy given by Eq. (j5j) determines the 
behavior of the coherence peak in a typical RIXS spec- 
trum in the SC state. An additional coherent excitation 
showing a linear dispersion down to uj = for q — > in- 
dicates the presence of nodes in the SC order parameter. 



FIG. 4: (color online) Dynamical structure factor of Eq. (3) 
in the main text of the paper for an isotropic Fermi surface 
with parabolic dispersion as a function of transferred momen- 
tum q near (0,0) for s wave (a) and d wave (b) superconduc- 
tors. The highest energy dispersion is drawn (dotted line) for 
gapped (a) excitations in the conventional SC state, while in 
the case of an unconventional superconductor (b) energy dis- 
persions (dotted lines) refer to excitations with momenta on 
nodal (gapless) and anti-nodal points (gapped) of the order 
parameter. 



In what follows we discuss the effect of the momen- 
tum dependence of the SC gap magnitude on the RIXS 
spectrum in the case transferred momentum q/kp <C 1. 
In this regime the spectral intensity is rather insensitive 
with respect to points where the sign of the SC order 
parameter changes. In the normal state (Ak = 0) the 
DSF is non-zero only for particle-hole excitations close 
to the Fermi surface leading to a coherence peak in the 
spectrum which disperses according to the Fermi veloc- 
ity. In Fig. 3] the dynamical structure factor is shown for 
a conventional s wave (a) and for an unconventional d 
wave (b) superconductor, calculated using Eq. (3) in the 
main text of the paper for an isotropic Fermi surface as 
a function of the transferred momentum q <C &f close to 
(0,0). Clearly seen is the absence of spectral weight for 
hu < 2 A in the conventional case (a), whereas at higher 
energies the coherence peak follows the dispersion £k of 
the bare electrons (dotted line). Instead, for a <i-wave su- 
perconductor (b) two types of excitations contribute to 
the low energy spectrum, related to two different regimes 
of momenta k contributing to the sum of Eq. (3) in the 
main text of the paper. If the momentum k is close to 
the nodal lines, i.e., the Brillouin zone diagonals shown 
in Fig. 1 of the main text of the paper, where Ak ~ 0, 
the excitations are particle- hole like with an energy ex q. 

The second type of excitations is provided by momenta 
close to the anti-nodal points k^ where Ak A = A is the 
maximum gap value. Here, for small values of q the 



Dynamical structure factor for strongly correlated 
electron system 

The main aim of this part of the supplemental ma- 
terial is to show that our main result for the phase de- 
pendence of the dynamical structure factor holds also 
in the presence of strong Coulomb repulsion of the con- 
duction electrons. To simplify the presentation of our 
arguments let us start from the limit of infinite strong 
correlations U oo where double occupancy is forbid- 
den l2[ . As is well-known using a Schrieffer-Wolf trans- 
formation the Hubbard model can then be replaced with 
a t- J model where the superexchange J 0, cf. [3]. Such 
regime of very strong correlations may be still regarded 
as rather realistic in describing many basic properties of 
the strongly transition metal oxides such as cuprates 0- 
[6[ . The superconductor is well described by the following 
Hamiltonian 



H 
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+ h.c). (6) 



The operators d\ a = d] a (l 



7t 



and d i(J = d ia (l 



n^-cr) are Hubbard creation and annihilation operators 
which enter since doubly occupancies of local sites are 
strictly forbidden due to the presence of strong elec- 
tronic correlations. They obey unusual anticommuta- 
tor relations. For instance, [di -,dt (7 ]+ = 5ijV a (i), with 
V a (i) = 1 - ni- a . 

For simplicity let us consider the dynamical structure 
factor X s (q, uj) for spin which is defined by the following 
correlation function 



X s (q,w) 



Jo 



d*<5i q (t)5*) w e 



i(uj-\-ir])t 



(7) 



Here, the time dependence and the expectation value are 
formed with Hamiltonian ((6]). To calculate the expec- 
tation value and the dynamical behavior, we diagonal- 
ize the Hamiltonian using new approximate quasiparticle 
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operators 7k, a > which are related to the original corre- 
lated d electron operators via = i£ k 7kt — ^kT^ki and 

d\tl = ^k7k| + ^kxLkf In the case of a sufficiently large 
hole concentration the operator V a (i) can approximately 
be replaced by its expectation value D and the Bogoli- 
ubov quasiparticle operators fulfill the following rela- 
tions [7]: [H^lj = E^, where £ k = +D 2 A^ 

and D = 1 — n/2. Replacing all d^ operators by the 

quasiparticle operators 7^, the time dependence can 
easily be evaluated. Each of the remaining expectation 
values contains a product of four quasiparticle operators. 
A final factorization leads to the following expectation 
values with two quasiparticle operators 



£k_ 



m k 



2EI 



(7k, CT 7k,<r) = -O- (7k,<T7k,<r), 



(8) 



where D = (V a (i)} = 1 — n/2 and and mk are defined 

by nk = (^ kc r^ kcr ) an d m k = (^kcr^kcr)- They are eval- 
uated using the Gutzwiller approximation (cf. Ref. [1), 
as nk = (D — q) + qf(e\<i) and mk — D — nk, where 
g = (1 — n)/(l — n/2) and /(sk) is the Fermi function at 
T = 0. 

Finally, we obtain the spin and charge dynamical struc- 
ture factor in the strongly correlated case 



xC ' S(q ' w) = 4iV X 



(9) 



k \ E k + £ k+q - (co + ir,) ^ q 47 k+q ,,)(7k,t7 k , t ) 
^(k,q ) 

— Ek — ilk+q — 



£ k - £ k+q - (w + irj) 



-— -y (7 k+q)t 7k +q ,t) <7 k4 7k,i> 
^<7 k+q ,t7k+ q ,t>(7k,t7 k ,t> J , 



where the coherence factors A + (k, q) and A (k, q) are 
denned by 



^ V ^k^k+q y 

(10) 



Hence, we conclude that the phase dependence of the 
dynamical structure factor in a strongly correlated elec- 
tron system is governed by the coherence factors given 
by Eq. (fT0]h which have, up to a constant D, the same 
form as in the uncorrelated case. 
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